Two-dimensional Vortices in Superconductors 
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Superconductors have two key characteristics. They expel magnetic field and they 
conduct electrical current with zero resistance. However, both properties are compro- 
mised in high magnetic fields which can penetrate the material and create a mixed 
state of quantized vortices. The vortices move in response to an electrical current dis- 
sipating energy which destroys the zero resistance state 1 . One of the central problems 
for applications of high temperature superconductivity is the stabilization of vortices 
to ensure zero electrical resistance. We find that vortices in the anisotropic supercon- 
ductor Bi2Sr2CaCu20g + a (Bi-2212) have a phase transition from a liquid state, which is 
inherently unstable, to a two-dimensional vortex solid. We show that at high field the 
transition temperature is independent of magnetic field, as was predicted theoretically 
for the melting of an ideal two-dimensional vortex lattice 2 ' 3 . Our results indicate that 
the stable solid phase can be reached at any field as may be necessary for applications 
involving superconducting magnets^ii^. The vortex solid is disordered, as suggested 
by previous studies at lower fields 7,8 . But its evolution with increasing magnetic field 
displays unexpected threshold behavior that needs further investigation. 



Shortly after the discovery of cuprate superconductiv- 
ity it was recognized that their high transition tempera- 
ture would mean that thermal fluctuations can produce 
a liquid vortex stated. In fact the thermodynamic tran- 
sition to superconductivity in a magnetic field occurs be- 
tween a thermally fluctuating liquid vortex phase to one 
or more solid phases™. The liquid vortex state is inher- 
ently unstable with non-zero electrical resistance. For 
extremely anisotropic materials, like Bi-22f2, the liquid 
phase covers a wide range of temperature; but it is not 
known exactly how wide this is in high magnetic fields. 
For fields perpendicular to the conducting planes, if He- 
axis, the transition temperature between liquid and solid 
vortex phases, T m (H), is principally controlled by vortex- 
vortex interactions which get stronger as the density of 
vortices increases, proportional to the field. The super- 
currents that form each vortex are mainly confined to 
the conducting planes and, in high field, they lose their 
coherence from one plane to the next so that vortices are 
expected to become two-dimensional. It is remarkable 
that the theory^ for the two-dimensional vortex melt- 
ing transition has only one parameter, the magnetic field 



penetration depth, and that this simple picture has not 
yet been experimentally confirmed. 

We use nuclear magnetic resonance (NMR) of 17 to 
detect the melting transition of vortices as a function 
of temperature and magnetic field. NMR can be per- 
formed on selected elements in site specific locations in 
the structur o n i 12 as we see in the 17 spectra in Fig. I. 
There are three stoichiometric positions for oxygen in Bi- 
2212, in atomic planes containing Cu, the 0(1) site; Sr 
planes, the 0(2) site; and Bi planes. The last oxygen 
atoms are unobserved^ owing to disorder in this plane. 
In addition there is a small amount of non-stoichiometric 
oxygen, 5-oxygen, too small to be observed directly by 
NMR and whose location in the structure is not yet 
established^. The parts of the NMR spectra in the mid- 
dle range of field in Fig. 1 are the central transitions from 
the 0(1) (broad) and 0(2) (narrow) sites. The electronic 
coupling to 17 is much stronger for 0(1) compared with 
0(2) as is apparent from the temperature dependences in 
Fig. 1. This is confirmed by measurements of the spin- 
lattice relaxation rate which are one order of magnitude 
larger for 0(1) as compared to 0(2). The narrowness of 
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FIG. 1: Spectra of 17 NMR in Bi-2212 for magnetic field 
parallel to the c-axis at fixed frequency. The sample is an 
overdoped, 28 mg, single crystal of Bi2Sr2CaCu2 0g+5, (Bi- 
2212), with T c = 75 K, with « 60% of 16 exchanged for 17 0. 
The optimally doped T c = 93 K. In field-sweep experiments 
(T = 20, 40, 60, 100 K), a decreasing Knight shift moves 
the NMR spectrum to the left. At high temperatures two 
oxygen sites can be distinguished. The central transition, 
(— | <-> +7f) for oxygen in the copper-oxygen plane, 0(1), is 
the wide line near 13.49 T at 100 K. The central transition for 
the oxygen in the strontium-oxygen plane, 0(2), is the narrow 
line near the zero Knight shift position at 13.51 T. The other 
two peaks are quadrupolar satellites of the O(l) resonance. 
For T= 4 and 12 K we show the Fourier transform of the 
echo at fixed field, having a smaller spectral bandwidth so that 
satellites are not observed here. Our separate measurement of 
the satellite intensity at T/T m = 0.67 compared to the central 
transition is the same as that at T/T m — 1.33, confirming this 
picture. The spectra at different temperatures are normalized 
to have the same peak intensity of the O(l) central transition. 
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FIG. 2: Temperature dependence of the 0(1) linewidth. Tem- 
perature dependence of the linewidth of the lr O NMR cen- 
tral transition presented as the square root of the second mo- 
ment of the frequency spectrum, 17 Ais, divided by magnetic 
field H in the range 3 to 22 T. Fields above 8 T were ob- 
tained at the National High Magnetic Field Laboratory in 
Tallahassee, Florida. Measurements at 27 and 29 T are in- 
cluded in Fig. 3 and 4, but omitted here for clarity. On 
decreasing the temperature there is a common behavior inde- 
pendent of magnetic field, represented by the smooth curve 
and its continuation to low temperature (dotted curve), that 
can be understood in terms of defects in the copper oxygen 
plane induced by oxygen doping. This curve is given by, 
Av = Avo + KHD/T, where Avq is a background magnetic 
contribution to the linewidth possibly associated with regions 
of the sample where superconductivity is suppressed, K(T) is 
the measured Knight shift, H/T is the field-to-temperature 
ratio and D is a Curie constant. The key feature of the data 
is the systematic break from this curve that we identify with 
the transition from a liquid to a solid vortex state. 



the 0(2) resonance indicates a homogeneous electronic 
environment with negligible spin-shift (Knight-shift) and 
uniform electric field gradient. For fast pulse repetition, 
as is the case for our measurements in Fig.l, 0(2) is par- 
tially saturated and comprises less than 6% of the spec- 
trum at 60 K, decreasing with decreasing temperature 
to 2% at 40 K. We remove it numerically in this range 
of temperature and below this, it has negligible contri- 
bution to the spectrum. The 0(2) resonance serves as a 
useful marker for the zero spin-shift position of 17 for 
both 0(1) and 0(2) sites. In this report we focus on 
17 0(1) NMR as a probe of the local magnetic field in the 
Cu02 plane. 

Decreasing the temperature below the superconduct- 
ing transition temperature we find the 17 0(1) resonance 
peaks move to the left, i.e. to a higher field at a 
fixed NMR frequency, (lower frequency at fixed field) ap- 
proaching the zero spin-shift position. Simultaneously, 
the central NMR line narrows. The decrease of the 



Knight shift to zero in the superconducting state is a 
characteristic signature for spin-singlet pairing. The line 
broadening with decreasing temperature in the normal 
state, Fig. 1 and 2, can be associated with a Knight 
shift distribution introduced by the <5-oxygen. Similar be- 
havior has been observed^ for chemical impurities, like 
Ni, Zn, or Li, substituted for Cu in the Cu02 plane in 
YBCO, and which form a local moment. Their contri- 
bution to the 17 O NMR linewidth is given by a Curie 
law, proportional to the ratio of applied magnetic field 
to temperature. Additionally, we find that below T c this 
linewidth is proportional to the temperature dependent 
Knight shift, which we measure independently, thereby 
accounting for the decrease with decreasing temperature 
in Fig. 2 as shown by the smooth curve, assuming a tem- 
perature independent residual contribution of 1.2 kHz/T. 
In this region, liquid vortex dynamics effectively average 
to zero their associated local magnetic fields. 
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On further cooling in the superconducting state there 
is a sharp onset for a new contribution to the linewidth 
which is not proportional to applied magnetic field. The 
well-defined temperature at which this line broadening 
appears is plotted in Fig. 3, decreasing, but progres- 
sively more slowly, with increasing magnetic field in our 
range between 3 and 29 Tesla. We identify this behav- 
ior with the formation of a solid vortex structure. More 
precisely, the extra linewidth corresponds to an inhomo- 
geneous magnetic field distribution that is static on the 
NMR time scale, » 0.1 ms and is asymmetric; see the 
T = 4 K spectrum in Fig. 1. This behavior is character- 
istic of the transition from liquid to solid vortex matter 
such as has been observed by /xSR— and SANS^ in Bi- 
2212, and NMR in YBCO 15 . Our observations are quali- 
tatively consistent with extension to high magnetic fields 
of the vortex melting phase diagram, H < 0.1 T, deter- 
mined by Hall prob o 16 i 17 and transport measurements^. 
Khaykovich et al*£ explore this behavior as a function 
of oxygen doping (anisotropy). NMR is complementary 
to these methods with definite advantages for detecting 
vortex melting at very high magnetic fields. The NMR or 
/iSR spectrum is a direct map of the magnetic field dis- 
tribution from vortex supercurrents. Specifically, NMR 
with 17 affords excellent resolution as a magnetometer 
on the atomic scale and has been exploited in previous 
work>i£ to spatially resolve and study excitations in the 
vortex core. 
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FIG. 3: The magnetic field-temperature phase diagram for 
vortex melting in Bi-2212 for H\\c. The transition tempera- 
tures, T m (H), are identified from the data, as shown in Fig. 
2. The two-dimensional vortex melting transition tempera- 
ture, T^f (vertical dashed line), is determined from a fit to 
the data. Deviations of the fit from the data are expected^ 
for H ~ H cr 



The strong upward curvature of the phase dia- 
gram in Fig 3. has been anticipated theoretically 3,10 . 
Torque magnetometry in fields H < 5 T gave similar 
indications^. For highly anisotropic superconductors the 



electromagnetic interaction between vortices dominates 
the Josephson coupling between planes. In this high field 
limit the simplest picture for vortex melting is a first or- 
der thermodynamic transition 3 - given by 



T m (H) = Tj 
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The limiting two-dimensional melting temperature in Eq. 
1. is, 



Tl D = A 



bid 



(3) 



where 0.605 < A < 0.615, comes from numerical cal- 
culations (Koshelev, A.E., private communication). This 
spread in A reflects the existence of an intermediate phase 
since the ideal two-dimensional melting scenario proceeds 
in two steps with an intervening hexatic phase bounded 
by continuous transitions 1 -. The melting temperature, 
Eq. 3, depends on the layer spacing 21 d = 1.5 nm, the 
flux quantum cf>o, Boltzmann's constant fee, and the sin- 
gle superconductive parameter A a h, which is the pene- 
tration depth for supercurrents in the Cu02 plane. The 
crossover field, H cr , depends mainly on the product of 
d and the mass anisotropy factor, 7 = \ c /\ a b — £,ab/£,c- 
Out-of plane components are denoted by a subscript c 
and £ is the coherence length. The numerical constants in 
the theory, Eq. 1. are b ss 1 and the exponent v = 0.37. 

The NMR vortex lineshape is asymmetric but less so 
than for a perfect line-vortex lattice suggesting that the 
vortex structure is somewhat disordered. In fact, /zSR 
and SANS show that at low temperature and H > 0.1 T, 
it is a vortex glass. We will assume that the difference 
in energy between disordered and perfect vortex struc- 
tures, at least in high magnetic field, is small compared 
to the energy for condensation from liquid to solid, and 
we use the framework of Eq. 1-3. to analyze the freezing 
that we have observed. From a fit to the data, dotted 
curve in Fig. 3, we find the cross-over field H cr = 2.5 T 
and 7 rs 78 and the two-dimensional melting transition 
temperature = 12 ± 1 K for which the penetration 
depth is X a i, = 220 ± 10 nm. Precision measurements of 
the absolute value of the penetration depth at low tem- 
perature are notoriously difficult. From earlier reports 
for Bi-2212, A a & is w 210 nm from magnetization data2i, 
269 ± 15 nm from cavity resonance methods 2 ^ and « 180 
nm from /1SR 7 . These results are in agreement with what 
we report here. From torque measurements, lye et al£^ 
found 7 > 200 for nearly optimum doped samples, and 
Watauchi et al£^ obtained 91 on overdoped material us- 
ing resistivity methods. Although 7 « 78 is in this range 
it must be considered only approximate since the theory 
is imprecise in the crossover field region. If we constrain 
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FIG. 4: Magnetic field dependence of the vortex linewidth and 
asymmetry. The vortex contribution to the second moment of 
the NMR spectra, ct(H), at low temperature, T ~ 4 K, is given 
in the lower panel as the difference between the data shown 
in Fig. 3 and the background contribution, smooth curves 
in Fig. 2, expressed as a second moment. A comparison 
with that expected from rigid-line vortices is calculated from 
Ginzburg-Landau theory 3 ^, smooth curve, taking the upper 
critical field to be 230 T. Low field /iSR measurement a 26 ' 27 up 
to f .5 T are shown as triangles with the lowest field less than 
H* in the Bragg glass phase, consistent with three dimensional 
ordering there, but dropping dramatically in the vortex glass 
phase. The linewidth data are normalized to the theoretical 
value at low field. In the upper panel we show the skewness, 
a = ([v — vi] 3 ) 1 /{[v — iA] 2 ) 1 / 2 where v\ is the first moment 
of the spectrum and the average is over the spectrum. A 
decrease in the second moment H w 6 T with no drop in 
asymmetry suggests intra-planar vortex ordering. The error 
bars are statistical, calculated from the NMR spectra. 



At low magnetic fields, H « H cr , it is well- 
established that there is a complex phase diagram for 
Bi-2212 with transitions on cooling from liquid vortex 
matter to a solid. At low temperature an increasing mag- 
netic field drives a transition at H* w 0.1 T from Bragg 
glass 25 to vortex glass, indicated by an abrupt increase 
in the symmetry of the /iSR spectrum 7 -, a decrease in 
its second moment 7 -, and disappearance of Bragg peaks 
in SANS 8 . At this transition vortices lose coherence be- 
tween planes. The resulting destructive interference be- 
tween vortices on adjacent planes averages out the mag- 
netic field distributio n 26 ' 27 ! 28 that one might expect for 
straightline vortices such as are seen in the Bragg glass 
phase 7 . This interference reduces both asymmetry and 
linewidth in the fj.SK spectrum. With even larger fields 
the second moment of the 17 NMR spectrum, cr(H), 
(ss linewidth squared) progressively increases and asym- 
metry is restored, until the spectrum linewidth collapses 
at 5 T. There are two possibilities for this anomalous 
behavior at 5 T. Either the transition is field-induced- 
ordering (ordering in-plane) or field-induced-disordering 
(further disorder between planes). The significant asym- 
metry in the spectrum in the high field phase compared 
to that at our lowest fields favors vortex ordering with 
increasing field. This scenario might follow if interplanar 
coupling is weakened with increasing field compared to 
intraplanar interactions, reducing frustration that origi- 
nated from vortices in adjacent planes. As a consequence 
vortices order two-dimensionally and we observe the sec- 
ond moment decrease abruptly. At much larger fields, we 
can see that the two-dimensional vortex solid becomes 
progressively more disordered. On theoretical grounds 29 
arbitrarily small amounts of impurity will lead to disorder 
of the vortex lattice. Sensitivity to quenched-in disorder 
should increase progressively with increasing magnetic 
field as the lattice spacing decreases with fixed vortex 
size. 



7, varying it over this wide range, our fit for is unaf- 
fected. For a less overdoped crystal, T c = 85 K, we find 
a slightly higher value, T~ = 16 K. 

Experiments consistent with the two-dimensional melt- 
ing theory have basic significance. But they are also rel- 
evant to applications of superconductivity at very high 
magnetic field. One of the most promising materials for 
magnet wire^^ to make a 30 T, superconducting, NMR 
magnet 6 is Bi-2212. To achieve this goal there will be an 
upper limit on the operational temperature of the mag- 
net determined by critical currents that are more easily 
stabilized in the vortex solid phase. 



To make a comparison with straight line vortices we 
have performed a Ginzburg-Landau calculation's, of the 
corresponding magnetic field distribution, Fig. 4. Our 
measurement is of a similar magnitude as expected from 
this calculation, although the GL approximation does not 
capture the details of the observed field dependence. In 
fact the transition from liquid to a solid in high mag- 
netic fields is predicted to occur first to a supersolid 
phase and then at lower temperatures to a decoupled 
solid phase where defects in the vortex lattice become 
zero-dimensiona l 10 ' 31 . NMR techniques may be helpful 
in exploring these fascinating aspects of vortex behavior 
in strongly anisotropic superconductors. 
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